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By constructing an Ising analogue of compact-star matter at sub-saturation density we explored 
the effect of Coulomb frustration on the nuclear liquid-gas phase transition. Our conclusions is 
twofold. First, the range of temperatures where inhomogeneous phases form expands with increasing 
Coulomb-field strength. Second, within the approximation of uniform electron distribution, the 
limiting point upon which the phase-coexistence region ends does not exhibit any critical behaviour. 
Possible astrophysics consequences and thermodynamical connections are discussed. 
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In absence of Coulombic interactions, the aspect of 
warm nuclear matter below its saturation density is a 
dense phase immersed in a low-density gas. The Coulomb 
repulsion strongly affects this liquid-gas phase transition 
by forbidding the condensation in macroscopic drops. 
This phenomenon, originally described in condensed- 
matter studies with the name of "frustration" , denotes 
all systems where no states exist where all interaction 
energies are simultaneously minimised: they range from 
magnets on specific lattices to liquid crystals, from spin 
glasses to protein folding. 

Nuclear matter in the density range of the liquid-gas 
phase transition can be produced either in the expansion 
and disassembling of nuclei involved in violent ion colli- 
sions, or in the core of supernovae explosions, as well as 
in inner neutron-star crusts. It is known since the seven- 
ties [l[ that the clusterised solid configurations in the in- 
ner crust of a neutron star give rise to complex phases , 
which are often quoted as pasta phases because of their 
suggestive topologies. Recent molecular-dynamics simu- 
lations show that these structures may survive also at fi- 
nite temperature ^ . The fluctuations connected to pasta 
phases are expected to enhance matter opacity to neu- 
trino scattering with important consequences on the su- 
pernova explosion and cooling dynamics Such a 
coherent neutrino-matter scattering is not only expected 
at low temperature, but even more in the possible occur- 
rence of a critical point in the post-bounce supernova 
explosion, with the associated phenomenon of critical 
opalescence [1, Q . At variance with usual- matter prop- 
erties, the expected increase in the static form factor was 
not observed in molecular-dynamics simulation of stellar 
matter at finite temperature @. This might be an ef- 
fect of the Coulomb interaction which also acts in finite 
nuclei. However, it should be observed that, differently 
from nuclei, star matter includes also electrons with the 
role of neutralising the net charge over macroscopic por- 
tions. 



The motivation of the present work is to provide gen- 
eral insights about the phase-transition phenomenology 
of neutral systems with Coulombic interactions. We 
construct an Ising analogue of compact-star matter by 
adding long-range interactions and a background of elec- 
trons to a Lattice-Gas model. Such an approach al- 
lows to take into account charge-density fluctuations by 
exact calculations, and to compare stellar matter with 
respect to other physical systems subjected to frustra- 
tion, like hot atomic nuclei^, @] and Coulomb-frustrated 
ferromagnets 9] . 

Let us consider a neutral system composed of charged 
particles occupying a cubic lattice of V cells with occu- 
pation numbers rii = 0, 1 immersed in a uniform back- 
ground of negative charge representing the incompress- 
ible degenerated gas of electrons 
charge qi = rii — n for each site i, where n 



IJ. This leads to a 



comes from the uniform background. 
Hamiltonian TYn+c = + 'He reads 
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where extends over closest neighbours. The isospin 
degree of freedom is not explicitly accounted, to allow 
a delocalization of the charge over the lattice The 
nuclear symmetry energy, here neglected, is known to 
change only quantitatively the phase diagram[llj. The 
short-range (nuclear- like T^n) and long-range (Coulomb- 
like Ti-c) interactions are characterised by the coupling 
constants e and Ae respectively, so that A measures the 
strength of frustration. To mimic stellar matter the 
numerical values are set to e = — 5.5MeT^ and Ae = 
ahcpQ^^x^ , where po — 0.17/m~^ is the nuclear satu- 
ration density, and x = 1/3 is a typical proton frac- 
tion. Tic can be rewritten as He = ^ Sij '^i'^jC'ij, 
where Cy = '"n^ ^ ^^^^ 
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= y J2}^o ^oi distance ry is imposed to be 

the shortest between i and j in the periodic space. To 
numerically accelerate thermodynamic convergence, the 
finite lattice is repeated in all three directions of space a 
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large number R of times, analogous to the Ewald sum- 
mation technique. Each site i has then R replicas of it- 
self, each one displaced from i of a vector mi, where m 
has integer components and L = \/V is the cubic lattice 
length. This procedure is equivalent to a renormalisation 
of the long-range coupling Cy . 

In ref. [9| a complete thermodynamic study of frustra- 
tion is dedicated to Ising ferromagnets, described by the 
Hamiltonian 



^Fc = 77 



(2) 



where Si = ni — 1/2 is a spin variable. The two models 
are related by: 

Wfc = Wn+c + 7T^A42 + - (3) 
2r 4 

where M. — ■'i ^^e magnetisation. We deduce 
that, contrarily to the standard form of closest-neighbour 
interactions, the charged Lattice-Gas model, eq.(IT|), in 
the grand canonical ensemble cannot be mapped into an 
Ising- ferromagnet model, eq.([2]), in an external magnetic 
field because of the term in Al^ in cq.([2]). Since this term 
scales with , a thermodynamic limit exists for Tipo only 
for vanishing magnetisation conversely to TYn+c |13l |. 
Since M., or the particle density p = A4/V + 1/2, is 
an order parameter of the Ising model, to impose = 
as a strict constraint can substantially modify the ther- 
modynamics of the system Moreover, the presence 
of the non- linear term directly affects the curvature 
of the order parameter distribution and so modifies the 
phase properties. 

Clarified these distinctions, henceforth wc analyse the 
thermodynamics of the TYn+c Hamiltonian. 

The multi-canonical ensemble [8] is a specifically 
adapted statistical tool to deal with frustrated systems. 
The Hamiltonian components TYn , Ti-c are treated as two 
independent observables associated to two Lagrange mul- 
tipliers (3rnPc, respectively. The multi-canonical parti- 
tion sum reads 

Zfi„MN,V) (4) 

where W{En, Ec,N, V) is the density of states with nu- 
clear energy E^, Coulomb energy Ec, and number of 
particles N. A generalized grand potential is defined by 



Zf3,,MN,V)e'''^'^''dN. (5) 



When f3c = (3n the ensemble coincides with the conven- 
tional grandcanonical form. When /3c = the system re- 
duces to the standard Ising model. Therefore, the multi- 
(grand)canonical ensemble allows to construct one single 



phase diagram both for neutral and charged matter by 
exploring the space (/3n,/3c)- It should be noticed that 
varying independently /?„ and Pc is equivalent to chang- 
ing the effective charge q^g = X(3c/ Pn- 

The canonical phase diagram in the (/3n, Pc) space can 
be accessed from the topological properties of the particle 
density distribution in the corresponding grand-canonical 
ensemble eq.© [l6|. For a given value of /3c, Fig- HJ^^ 
describes, as a function of the inverse temperature /3n; 
the grand-canonical density distribution calculated at the 
critical chemical potential ii — iic- This distribution is 
sampled from eq.® with a standard Metropolis tech- 
nique [l3|. At high Pn, the density distribution shows 
two peaks of the same height, which are associated to 
two coexisting phases [3 ■ These two phases join at the 
limiting temperature (point (4)) above which the distri- 
bution stops to be bimodal. From a series of similar cal- 
culations performed for different values of /3c, we could 
extract the evolution of the inverse limiting temperature 
as a function of /3c as plotted in Fig. [TJ3. We observe an 
increase of the limiting temperature /3~^ for increasing 
strength of the Coulomb field (about 6% for a proton 
fraction x — 1/3). 

It is important to notice that many other physical sys- 
tems subjected to Coulomb frustration exhibit the op- 
posite behaviour. This is notably the case of frustrated 
Ising ferromagnets, as well as of finite atomic nuclei. In 
such cases the Coulomb repulsion is known to reduce the 
limiting temperature 0, H, 0, [13, [111- This reduction 
is also an usual expectation in the astrophysical con- 
text [1, [3, [Ifl,] . However, a recent calculation of nuclear- 
pasta structure within the RMF model [I4I indicated a 
widening of the density range connected to the mixed 
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FIG. 1: (Color online) Metropolis calculations (/i — fic, L = 
10): (A) The logarithmic cluster plot gives the density distri- 
butions for different values of /3n and for /3c = 0.83e, which is 
the inverse limiting temperature computed for Pc = (3n- The 
solid line indicates the ridge of the distribution, the dashed 
line the coexistence region of an Ising-like system (/3c = 0). 
The points (2) and (4) indicate the limiting temperatures for 
Pc = and Pc = Pn- (B) Phase diagram in the multicanon- 
ical ensemble giving the inverse limiting temperature Pn as a 
function of Pc- The points (1) and (3) are used in Fig. [5] 
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phase region when the Coulomb field is included under 
the constraint of charge neutrality over the Wigner-Seitz 
cells. Our results are consistent with this finding, and 
suggest that this effect is model independent and should 
persist at finite temperature [l^. 

The analysis of the event distribution elucidates the 
physical origin of the increase of the limiting tempera- 
ture. The left side of Fig. [2] shows the transition from 
subcritical to critical partitions for an uncharged system 
(/3c — 0). When the system is uncharged, the observ- 
able Ec defined in eq.([T]) does not represent a physical 
Coulomb energy, and measures the compactness of the 
system. Its contribution does not influence the partition 
probability since /3c = 0, and compact clusterised parti- 
tions with high Ec can be explored. 

When the system is charged (right side of Fig. [2]), Ec 
enters in the partition distribution eq.® as the physical 
Coulomb energy. This energy is maximal in fragmented 
configurations at p/po ~ 0.5, and minimal in pure-phase 
events, where the largely uniform proton-charge distri- 
bution is almost exactly compensated by the electron 
background. The pure-phase events are then increasingly 
favored with increasing charge. Thus, bimodal distribu- 
tions are still found in the charged system (panel 3) at 
temperatures at which critical events dominate the un- 
charged system (panel 2). This is entirely due to the 




EJV[e] 



Contours: L = 10 
Clusters: L=20 



Er/V[e] 



FIG. 2: (Color online) Probability distribution for the two 
energy components En,Ec for the four points indicated in 
Fig. [l]3. Contour plots refer to L = 10; logarithmic cluster 
plots refer to L = 20. Points (1) and (3) belong to the co- 
existence region and manifest bimodal patterns. Both point 
(2) and (4) are very close to the limiting temperature but 
the critical scaling effects exhibited at point (2) are absent at 
point (4). 



screening effect of the electrons, while in finite nuclei 
the highest Coulomb energy is associated to the homoge- 
neous liquid-like partitions, leading to the opposite effect. 

The question now arises as to whether the limiting 
points of the frustrated system are second-order critical 
points like in the uncharged system, or in the Coulombic 
RPM model describing phase separation in electrolytes 
{211] . Criticality arises from the divergence of the corre- 
lation length ^. It describes the exponential decay of the 
correlation function (T(rij) = (SniSrij) = (ninj) — (ni)(nj) 
according to the expression a{r) cx e"*"/^ • r~^^~^+''\ 
where D is the space dimension and 77 is a critical expo- 
nent. Using the properties of the Cy matrix, a(r) can be 
easily related to the mean Coulomb-energy density by 



/Ec 
\ V 



aiO) 



Ae X ^ 



(6) 



Eq. ^ indicates that a diverging correlation length at the 
critical point manifests by a divergent Coulomb-energy 
density. When /3c 7^ the only solution for the system 
to avoid such a singularity is to suppress the second- 
order critical character of the limiting point [l^ . In this 
case, the correlation length keeps a finite value at the 
limiting temperature and the thermodynamic limit is ful- 
filled by a constant value of {Ec /V) . Such a behaviour is 
tested in Fig. [21 where the same calculation is repeated 
for different lattice sizes: the average Coulomb energy is 
seen to increase with the lattice size in the Ising system, 
while Coulomb-energy fluctuations appear quenched in 
the frustrated system. It is interesting to note that this 
argument does not apply in other neutral Coulombic sys- 
tems 2l|, where density fluctuations do not necessar- 
ily imply charge fluctuations, and can therefore diverge 
keeping a finite Coulomb energy. For the same to be 
true in the proto-neutron star, the electron field should 
be strongly polarized at the limiting temperature, lead- 
ing to a complete charge screening of the dishomogeneous 
pasta structures. Due to the high incompressibility of the 
degenerate electron gas, this is likely to be unphvsicalfl^. 

The quenching of criticality can be formally verified in 
terms of critical exponents. If the asymptotic value of 
the limiting temperature Tum = limi^^oo 7iini(i) corre- 
sponds to a critical point, finite-size scaling insures that 
Tiinu(i) should evolve as fiim{L) - Tum oc L-^/" 
Fig. [3] illustrates that, while the scaling is respected with 
the Ising value of v when /3c = 0, should be in- 
finitely large to describe the frustrated system. Since 
ly rules the divergence of the correlation length, ^ cx i"" 
with t = T/Tiim — 1, this is a first indication of a fi- 
nite correlation length for the charged system. We then 
test the finite-size scaling theory [2^ on the quantity 
X — '^ij/T ■ This form of x represents the suscep- 
tibility only for T > T\nn{L). In this case, it should scale 
with the critical exponent 7 as L^^'l^x = fiL^^^'^t) if 
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the limiting point represented a critical point; the func- 
tion f[L^^/^t) should be constant when ^ ~ L and 
/(L-i/'^t) cx (2.1/"*)-^ when ^ < L. When T < rum(L), 
X contains also jumps between the low-density, (n)G, 
and the high-density, (n)L, solutions and should therefore 
scale with the critical exponent /3 defined from the scaling 
of the order parameter {n)i^ — {n)Q oc as x 
By introducing the hyperscaling relation D = {^ + 2/3) /i/, 
all scaling laws condense in the form 



-7/1/ 



X oc 



constant 



as ^ ~ L 

as ^ « i, T > 

as ^ « i, T < Tiim(i) 



(7) 



Fig. [3p illustrates the perfect consistency of the /3c = 
system with Ising critical exponents. Conversely, for the 
charged system, no combination of /?, 7 and u can be 
found to let all calculated points collapse on the same law. 
In particular, for 7/1^ = and — > 00, only the points 
for T ^ Tiim(L) collapse on one line, while those for 
T <^ T\im{L) disperse as shown in Fig. [3jl). This demon- 
strates that the effect of the Coulomb field is not a simple 
increase of the v exponent, but a complete quenching of 
criticality for the frustrated system. The loss of critical 
behaviour has been already observed in Ising models for 
the Ising frustrated ferromagnet, eq.(l2]), where the coex- 
istence region was seen to end at a first-order point ^ . 
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FIG. 3: (Color online). Panels A,B. Evolution of Tum as 
a function of the linear size L, and extrapolation towards 
the thermodynamic limit. Panels C,D. Study of the finite- 
size scaling eg.©. All calculations are presented for the un- 
charged (/3c = 0) and the charged system (/3c = /3n) with 
proton fraction x = 1/3. 



Our conclusion is twofold. First, we observed that, in 
presence of a uniform electron background, the Coulomb 
field originating from charge-density fluctuations in- 
creases the limiting temperature for phase coexistence, 
at variance with the strong decrease obtained in the ab- 
sence of electrons or in mean-field calculations [lil [20j . 
Therefore, the mixed-phase phenomenology may be rel- 
evant for the proto-neutron-star structure up to slightly 
above the critical temperature of normal nuclear mat- 
ter (Tc w 15 MeV for symmetric matter), in a wider 
temperature-range than usually expected [2]. Second, 
the Coulomb field suppresses the critical character of the 
limiting temperature. For this reason we expect warm 
stellar matter to show small opacity to neutrino scatter- 
ing in agreement with the findings of ref . (23j . 
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